Bound-states in dimerized and frustrated Heisenberg Chains 
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Using the Bond Operator Technique (BOT), we have studied the low energy excitation spectrum 
of a frustrated dimerized antiferromagnetic Heisenberg chain. In particular, we have compared our 
analytical results with previous Exact Diagonalizations (ED) data. Qualitatively, the BOT results 
are in good agreement with the ED data. And even a very good quantitative agreement is obtained 
in some parameter region. It is clearly shown that there is only one elementary excitation branch 
(lowest triplet branch) and that the two other well defined excitations which appear below the 
continuum, one singlet and one triplet, are bound states of two elementary triplets. 
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I. INTRODUCTION 

The magnetic properties of low dimensional quantum 
spin systems have attracted a considerable attention of 
theorists over the last decades. The interest for these sys- 
tems was boosted by the recent discovery of non organic 
spin-Peierls materials like CuGeO^ 0] and A^aV205 [HJ|] 
and spin ladder compounds like SrCu20z, Sr2Cu30^ 
It turns out that some of the properties of these 
materials can be described by the frustrated dimerized 
Heisenberg model. In the one dimensional case, the 
Hamiltonian of this model reads. 



if = ([1 + (5(-l)'-]S, • S,:+i + aS, • S 



1+2 



(1) 



where i denotes the sites of a chain with length L and 
Si are S = 1/2 spin operators. J > is the nearest- 
neighbor exchange coupling, a the frustration parameter 
from next-nearest neighbor coupling and 5 is the dimer- 
ization parameter. 

Depending on the coupling parameters, this model ex- 
hibits various phases including spin-liquid phase (gapless 
phase), dimerized phase, Neel ordered phase and incom- 
mensurate phase. The ground-state properties of this 
model were extensively studied using various numerical 
and analytical tools, like Exact Diagonalizations method. 
Quantum Monte Carlo, DMRG (density matrix renor- 
malization group), Bethe ansatz, bond operator tech- 
nique, etc.. Another interesting feature of the frustrated 
dimerized Heisenberg model is the possible existence of 
well defined excited states below the continuum (bound- 
states?) ^-|^. Recently, Bouzerar et al. have provided 
a detailed numerical analysis (Exact diagonalizations 
method) of the excitation spectrum in this model [||. 
It was shown that: (i)The Singlet- Triplet excitation gap 
ratio is a universal function which depends on the frustra- 
tion parameter only and (ii)A well defined second triplet 
branch appears below the continuum in the vicinity of 
momentum q = Q (resp. tt) only if the strength of the 
frustration is large enough. A useful approach to describe 



the dimerized phase is the Bond-Operator representation 
of spins introduced by Chubukov and Sachdev and 
Bhatt . Recently a number of works have been carried 
out with the help of bond-operator method. Employ- 
ing the bond-operator methodology Brenig studied 
the effect of interchain coupling in spin-Peierls systems 
like CuGeO^. This method was also successfully used by 
Sushkov and Kotov to spin ladder model in the strong in- 
terchain coupling limit » J//, to show the existence 
of a singlet and triplet bound states below the contin- 
uum [|l4| . In this paper we will study the excited states 
of the frustrated dimerized Heisenberg model using the 
Bond-Operator Technique (BOT) and compare the re- 
sults of our calculation with the Exact Diagonalizations 
data (ED) ||^. We will show that the agreement between 
ED and BOT is very good at least qualitatively, and even 
quantitatively in the vicinity of momentum q — 7r/2. We 
will also show that the lowest triplet branch is the only 
elementary excitation and that the singlet and the triplet 
which appear below the continuum are bound-states of 
two triplets. 



II. BOND OPERATOR TECHNIQUE 

In this method the Hilbert space of the spin degrees 
of freedom is represented in terms of singlet and three 
triplet states which is created out of vacuum by one sin- 
glet [s^) and three triplets (i^) creation operators 
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A representation of spins in terms of the singlet and 
triplet operators is given by 
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where a, (3 and j = x,y, z and eaf3-y is the Levicivita 
symbol representing the totally antisymmetric tensor. 
Henceforth it is assumed that all repeated indices over 
a, f3 and 7 are summed over. The four operators (s, t^, ty 
and tz) satisfy the usual bosonic commutation relations. 
In order to ensure that the physical states are either sin- 
glet or triplets one has to impose the constraint 



1. 



(4) 



It is a difficult task to fulfill this constraint exactly (in 
each bond) and in most of the works this constraint is 
satisfied only in average. To overcome this difficulty, an 
infinite on site repulsion has been introduced by Kotov 
et al. [|3|. In Bond-operator representation the Hamil- 
tonian reads, 

H = Ho + Hi+H2 (5) 



9 ^ f 

H = --J{l + d)L+ J2 ^'^(aLa^fe+o) (6) 
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with Afe = J(f + (5) - 4(1 - 5 - 2a)cos{2k) and Bk = 
— 1(1 - 5 ~2a)cos{2k). 

Let us now take into account the effect of the local 
constraint eq. (H). Following the ref. (l| we include 
in the previous Hamiltonian an infinite on site repulsion 



(8) 



= -| J(l + 5)^ + J(l + <5) E 



- ^0 ^{ta^to,^+l4si+l + ^L^L+l + h.c) 

i=l 
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with Kq 
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2a),K, ^ {{1 



2a) and 



K2 = i{l-d). 

This method gives the exact ground-state in the limit 
when the ground state wave function consists of a prod- 
uct of local independent dimers which is especially real- 
ized on the 'disordered line" 2a + S = 1. On this line the 
Ground-state energy per site is Eq/L = — |(1 + 5). It 
should be noted that in the BOT the elementary excita- 
tion is a local triplet (nearest neighbor sites). 

In order to perform the calculation of the triplet dis- 
persion: (i) Initially we neglect completely the local con- 
straint which is equivalent to set Si = 1. (ii) We restrict 
ourself to the quadratic terms only. As it was previously 
shown that the effect of the higher order terms Hi and 
H2 are small [|5[0]. 

After using the well known Bogoliubov transformation 
tak = 'U'kO'ak + '^k^'a-k^ ^^"-^ hamiltonian reduces to, 



Concerning the calculation of the self-energy correction 
to the dispersion we follow the diagrammatic approach 
developed in ref. p^ . Let us just summarize the main 
steps of the calculation. At first the vertex scattering am- 
plitude r(A:, oj) is evaluated by solving the Bethe-Salpeter 
equation shown in fig.0 a. 

r{k,uj) = [i / dqduj'G{q,Lj')G{k — q,ui ~ uj')]-^ 



N UJ - UJ„ - LUk-q 

q ^ ^ 



+ {u^v,uj^ -u}]-^ (9) 



where G{k,uj) is the normal Greens function of Hq. 
We neglect the higher order terms of Vq since ^qV^ is 
proportional to the number of triplets in the system and 
it is very small in the dimerized phase. 

Then taking the scattering vertex as T{k, oj), the corre- 
sponding self energy is obtained from the diagram shown 
in fig I b. 



^(.k,uj)^— vlr{k + q,uj - ujq). 



(10) 



1 = -2 



Now we calculate the renormalized spectrum 

nk = Zfc^((Afc+S(fc,0))2-B2) (11) 

from the poles of the modified Greens function. Here the 
quasiparticle weight 



Zk 



(12) 
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III. SINGLET-TRIPLET GAP AND DISPERSION 

In this section we will compare the BOT results (with 
and without the constraint) of the Singlet- Triplet gap and 
Dispersion to the ED data of ref . Q . We observe in fig. 
that the gap A bot is always larger than the exact result 
and even gets worse when the constraint is included. For 
a given 5, the difference between the gaps with and with- 
out constraint reduces when it approaches the disordered 
line. On the disordered line, it can be seen easily from 
equation ( [ll| ) , that the dispersion is unchanged when the 
constraint is included, since E(fc,Li;) cx (2q! + (5— 1). How- 
ever, the agreement between BOT and ED gets better 
in the limit of large 5. For instance at a = 0, we get 
-fe^ = 0.67 for 5 = 0.2, but this ratio is 0.87 when 
— 0.4. This ratio increases with a. Furthermore, the 
curve including the constraint is parallel to the ED data 
curve, this means that the effect of frustration is taken 
into account more properly when the constraint is in- 
cluded. However, in the limit of vanishing 5 the effect of 
frustration is not taken properly in the BOT approach 
For instance at the Majumdar-Ghosh [|l6| point 
(a ^ 0.5,(5 = 0.), -fe^ « 0.25. 

Let us now discuss the dispersion of the triplet ex- 
citation {u){q)). In order to analyze the effect of the 
frustration on Lo{q) we have fixed the dimerization pa- 
rameter 5 = 0.2. In fig. ^ a (resp|| b) we have plotted 
the dispersion for a = (resp. a = 0.3). In fig|| a 
and figjl b we observe that the agreement between the 
ED data and BOT with constraint are excellent in the 
vicinity of q = 7r/2. This agreement will get better with 
increasing 5. Note that on the disordered line the BOT 
data coincides exactly with the ED data. Indeed, on 
this line and at q = 7r/2, the lowest triplet excitation 
is |T(7r/2) >= Y.I e'^'-Sl\GS >, where the ground-state 
wave function \GS > consists of product of independent 
dimers [|l^. Note that even far from the disordered line, 
this state remains a good approximation of the exact low- 
est triplet state in the vicinity of tt/2. It is important to 
emphasize, that away from q = n/2, the disagreement 
between BOT and ED data is mainly due to the fact 
that the triplet excitation are not local objects, but in 
the BOT scheme by construction the triplet excitation 
are local. We observe in fig.|| a (resp.|| b) that when ap- 
proaching the disordered line the width of the dispersion 
strongly decreases, and the agreement with the ED data 
gets worse. On the disordered line, the BOT triplet ex- 
citation is dispersionless and uj{q) = J{1 + S). However, 
as discussed previously, when S is increased the agree- 
ment between BOT and ED gets better and better and 
coincide for the special case 6 = 1 and a = 0. 

IV. TWO MAGNONS BOUND-STATE. 

The quartic term Hi in the Hamiltonian eq. ^, con- 
sist of triplet-triplet interaction. This term can lead to 



an attractive interaction between triplets to form bound- 
states. This bound-states can be singlet, triplet or quin- 
tuplet. We will focus only on the possibility of a singlet or 
triplet bound-states which were observed below the con- 
tinuum in the ED calculations Q . In the previous section 
we have seen that the effect of the constraint does not sig- 
nificantly modify the triplet dispersion, thus to simplify 
the calculations we will neglect the constraint. 

A. Singlet bound-state. 

To perform the calculation of the singlet bound-state 
, let us consider the most general singlet wave function 

l^'iQ) E *''(fc)«lQ/2-fc«lQ/2+fe|0 > (13) 

k,a 

of total momentum Q, constructed from two triplet ele- 
mentary excitations, where $"^(/c) is determined from the 
Schrodinger equation, 

i/|*'^(Q) >=£;^(Q)|*^(Q) > (14) 

with singlet bound state energy E^{Q). This leads im- 
mediately to the integral equation of the form, 

= -4KMq) r -9{P)<^^{P) (15) 

where g(jp) = cos{2p)uQ/2-pUQ/2+p- 

In the vicinity of the disordered line (|2a-|-(5— 1| << 1) 
this integral can be solved analytically. In this region 
of parameter Ufe 1 and uik/ J ~ (1 + (5) — ^(1 — (5 — 
2a)cos{2k). The corresponding singlet energy for the 
momentum Q is, 

E^{Q)/J = 2{\ + 5) 

+ 2a)2 + 4(1 - ^ - 2a)2cos2(Q). (16) 

For any momentum Q the energy of the singlet state 
is always smaller than the energy of the lower edge 
of the continuum defined as E^{Q) = minq{ujQ/2-q + 
ujQ^2+q) = 2(1 -h <5) - 1(1 - (5 - 2a)cos{Q)\. Thus, a well 
defined singlet bound-state exists for any momentum Q. 
Far from the disordered line the equation (|l^) should be 
solved numerically. 

B. Triplet bound-state. 

Let us now perform analogous calculations in the 
triplet sector. The most general triplet wave function 
can be written, 

l*a(Q) >- E *«(fc)^«/374Q/2-fc<Q/2+JO > . (17) 
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The corresponding Schrodinger equation is, 



2Kihiq) 



dp 



hip)^l(p) (18) 



where h(p) = sin(2p)uQ/2„pUQ/2+p. 

As we did it previously, we solve this equation analyt- 
ically in the vicinity of the disordered line. We get, 

E^iQ)/J = 2(1 + d)-K,[l + i'^fcos'iQ)] (19) 

Ai 

We find that the Triplet is below the continuum (i.e. 
E^{Q) < E^iQ)) when Q > Qc, where 



_i l-5 + 2a 
Qc — cos [- 



'2(1 -S -2a)' 



(20) 



C. Discussions. 

In this subsection we will compare the BOT results to 
the ED data. Let us first discuss the dependence of the 
singlet-triplet ratio R — where A'^ is the singlet- 
singlet gap and the singlet-triplet gap (A-^ = A). It 
was shown in ref. |^ that, in presence of dimerization 
((5 0) this ratio is a universal function which depends 
on the frustration parameter only (in the commensurate 
region). In fig. ^ we have plotted R for different param- 
eter (5 as a function of a. The agreement with the ED 
data is surprisingly good. We observe that R^^^ has a 
small dependence on 5. However, when 5 is increasing 
j^BOT _^ R^^ , Especially for a = we observe that 
R{a — Q)^2. The deviation for large 5 in the ED data 
were attributed to the crossing of the disordered line. 

It was shown in ref. that a well defined second 
triplet branch split from the continuum in the vicinity 
of g = (resp. q = tt) li the strength of the frustration 
is large enough, a > a* {5). In particular it was shown 
that a* {5) « (1 - (5)/3. In fig.^ we have plotted a* as 
a function of 5 calculated exactly and within the BOT 
method. We find that the qualitative agreement with the 
ED data is very good. For large 5 we have found that 
o^bot(^) ~ l/2a^£i(i5). We believe that this discrepancy 
on the slope is due to the fact that the width of the 
dispersion of the elementary triplet is underestimated in 
the BOT approach. 

We observe that in the unfrustrated case {a = 0) in- 
creasing 5 reduces the region where the second triplet ap- 
pears below the continuum. In other words, Qc{5), the 
momentum where the triplet split from the continuum, 
increases with i5. The variation of Qc with 5 is plotted in 
figj^ and it is in agreement with eq.(20), Qc — > 7r/3 when 
5^1. Thus, in the unfrustrated case, the effect of 5 
reduces the effective interaction between two elementary 
triplets excitation. In order to visualize the effect of the 



frustration, we have plotted in fig.0, the dispersion of the 
lowest excitations, ie. the elementary triplet excitation 
and the two bound-states for a fixed value of the dimer- 
ization parameter. The figure shows, in agreement with 
the previous work ||], that the effect of the frustration 
increases the region where the second triplet is a well de- 
fined excitation. These figures are very similar to fig. 5 of 
ref. §. 

V. CONCLUSIONS. 

As a conclusion, using the Bond operator method we 
have shown that in the frustrated dimerized Heisenberg 
model, there is only one elementary excitation branch 
(lowest triplet branch). Depending on the parameters 
{5, a), two bound states of these elementary excitations, 
one singlet and one triplet, can appear below the contin- 
uum. Even in absence of frustration the singlet bound- 
state is a well defined excitation for any momentum. 
However, in good agreement with previous work, the re- 
gion where the second triplet is well defined depends on 
the frustration parameter. The triplet bound-state is ob- 
servable in the vicinity of g = (resp. q = tt) when the 
frustration strength is large enough. Furthermore, in the 
unfrustrated case, the region where the second triplet is 
split from the continuum, reduces when increasing the 
dimerization parameter 6. 
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FIG. 1. Diagrams in the ladder approximation for the ver- 
tex scattering amplitude (a), and the self-energy (b) 



FIG. 2. Gap A as a function of a for 5 = 0.2 and 0.4 
.The full symbols are Exact Diagonalizations data from ref. 
(extrapolated in the thermodynamic limit). The continuous 
line are BOT calculations without constraint, and the dotted 
line with constraint. 



FIG. 3. Dispersion of the lowest triplet branch for a fixed 
5 = 0.2 and a = (a) and a = 0.3 (b). The symbols are ED 
data (for a system size L = 20). The dashed line are BOT 
calculations without constraint, and the continuous line with 
constraint. 



FIG. 4. Singlet-Triplet gaps ratio. The symbols are ED 
data from ref. The lines are BOT data calculated for 
different values of S. 



FIG. 5. a* as a function of S. Symbols are ED data from 
ref. M, and the contiimous line corresponds to BOT results. 



FIG. 6. Qc as a function of 5 in the unfrustrated case. 



FIG. 7. Dispersion of the triplet elementary excitation, sin- 
glet bound-state and triplet-bound-state for a fixed 5 — 0.2 
and Q = 0.1 (a), 0.17 (b) and 0.25 (c). The continuum corre- 
sponds to the shaded region. 
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